Atom-to-molecule conversion efficiency and adiabatic fidelity 
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The efficiency of converting two-species fermionic atoms into bosonic molecules is investigated in 
terms of mean-field Lagrangian density. We find that the STIRAP technique aided by Feshbach 
resonance is more effective than the bare Fechbach resonance for 6 Li atoms rather than 40 K atoms. 
We also make general consideration on the symmetry and its relevant conservation law, which enable 
us to introduce a natural definition of adiabatic fidelity for CPT state. The calculated values of 
the fidelity then provide an interpretation on why the conversion efficiencies for 40 K and 6 Li are 
distinctly different. 
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I. INTRODUCTION 

There has been much attention to the molecular Bose- 
Einstein condensate (BEC) which is a versatile play- 
ground not only for cold atomic physics experiments but 
also for other research areas, such as condensed matter 
theory as well as quantum field theory. It is no longer 
a pure Bose system as the molecular BEC can be a sys- 
tem of Bose-Fermi mixtures. This makes the theoreti- 
cal structure rich and colorful. To realize the molecular 
BEC, one must create stable molecules with long lifetime. 
In the recent experiments [2, H, IH , the technique of Fes- 
hbach resonance plays an important role in the creation 
of molecules. Since the molecules created through such a 
technique usually suffer from fast decay due to the vibra- 
tional excitation, the stimulated Raman adiabatic pas- 
sage (STIRAP) in photoassociation 0] has been regarded 
as an effective approach to create ground-state molecules. 
The success of STIRAP technique requires the existence 
of a coherent population trapping (CPT) state which can 
be followed adiabatically Such a condition can be 
fulfilled for linear A systems by appropriately choosing 
laser frequencies. However, for the system with inter- 
particle interactions, the two-photon resonance condition 
dynamically changes when population is transferred from 
atomic states to molecular states. This makes the CPT 
state more difficult to be followed adiabatically. The adi- 
abatic property was studied by means of adiabatic fidelity 
in a recent theoretical work Q for a simplified model 
of monoatomic system without inter-particle interaction. 
It is worthwhile to appropriately define the adiabatic fi- 
delity to study the adiabatic property for a more realistic 
system, such as the diatomic system which has been re- 
alized in serval experiments @,H, H, EL E2| • 

In this paper, we consider systems consisting of 
fermionic atoms in different hyperfine states and their 
compounded molecules coupled through the STIRAP 
technique aided by Feshbach resonance. We discuss sys- 
tems of 6 Li and 40 K as concrete examples. In section ILT1 
we model the systems with inter-particle interaction and 
derive the dynamical equations for them. We make a 
general consideration on the symmetry and the relevant 



conservation law, and then introduce the definition of fi- 
delity for our system. In section !!!!! we look for solutions 
of the CPT state. In section |PVl we calculate the conver- 
sion efficiency for 40 K atoms and 6 Li atoms, respectively, 
and discuss the corresponding features. In section [Vj we 
study the relation between the atom-to-molecule conver- 
sion efficiency and the adiabatic fidelity for CPT state. 
We also study the effect of the decay of quansibound 
molecules and compare the difference between K and 
Li systems. Our results are briefly summarized in sec- 
tion EH 



II. MODEL AND ITS GENERAL PROPERTIES 

We consider that two species of fermionic atoms are 
converted into stable molecules via the STIRAP tech- 
nique aided by Feshbach resonance. Here we use | a) and 
| b) stand for the states of the free atoms in the open 
channel, and | m) and | g) for the quasibound and ground 
molecular states in the close channel, respectively. A 
pair of ferimonic atoms are coupled with the quasibound 
molecular state | m) through a Feshbach resonance, where 
the coupling strength is denoted by a' and the detuning 
by £' . Additionally, the sates | m) and | g) are coupled 
with each other through a laser field with the coupling 
strength SI' and detuning A'. Since the system is a Fermi- 
Bose mixture whose Hilbert space actually carries out 
the representation of a graded unitary group SU(2 | 2), 
it will be an arduous work to study the dynamics of the 
system by means of the Heisenberg equation of motion. 
Whereas, with the help of mean- field Langrange den- 
sity fl3j , one can investigate the dynamics conveniently, 
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where T, denotes the kinetic energy term of the zth 
component and i,j — a,b,m,g. Here the coefficients 
A'., = A', = 2irh 2 an/m i4 are the interaction strengths 

tj J t ij / tj U 

between particles with — m^m ■/ (m i +m^J being the 
reduced mass and the s-wave scattering length. 



A. The conservation law 

The Lagrangian ([Ij is no more invariant under a si- 
multaneous global phase transformation due to the pres- 
ence of the atom-to-molecule conversion term. As the 
Lagrangian does not include the term flipping the two 
fermionic components into each other, their correspond- 
ing phase parameters are not necessarily the same in gen- 
eral. One can find that the Lagrangian Eq. ([TJ) is invari- 
ant under the phase transformation, (if) a , ip b , ip m , ip g ) i— > 

(■lpa,1pb,i> mi 'lpg)U('8a 1 'd b ) with 
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Following the formulism of Noether theorem, we evalu- 
ate the variation of the action caused by the infinitesimal 
phase transformation U(5d a , fi$ b ), 



51 = J {L(ip+S^,^* + Sp,8^+ 88^,8^* + 88^*) 

For simplicity in the above expression, we omitted the 
subscripts of ip labelling different components, abbrevi- 
ated d{x) for dtdxdydz and for d/dx^ with {x^} = 
{t, x,y, z}. Because of, explicitly, 5ip a — i(S'&a)4'a, 8tp b — 
i(dtf b )ijj b , Sipm^g) = i(^o + ^i?b)^ TO ( g ), and their complex 
conjugations, we derive the following result, 

SI = J {{d^Sda + (dMSA}d(x), (3) 

where J% = J£ + + J£, j£ = J% + J£ + J£ with 



J? = i>i 
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— , i = a, b, m, g. (4) 
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The system is invariant under the transformation Eq. @ 
such that SI = 0, which gives rise to two conserved cur- 
rents. 



9 M J£ = 0, 8^ = 0. 



(5) 



In the present paper, we focus on uniform system 
(Vtpi « 0) and hence neglect the kinetic and trapping 
potential terms. Then the conservation law Eq. {5]) reads 
d , io i io i , i o \ „ d 

dt' 

0, which implies that 



dt -(l^| 2 + IV'm| 2 + |^| 2 )=Q, ^(l^| 2 + l^n| 2 + |^| 2 ) 



m 2 + \^ m \ 2 +m 2 = n a , 

\A\ 2 + \^m\ 2 + \i>g\ 2 = n b , 



where the constants n a and n b are determined by V, the 
volume of the system, together with N a (0) and N b (0), the 
initial numbers of species a and 6, i.e., n a = N a (0)/V, 
n b = N b (Q)/V. Here we assume there are no molecules 
at the initial time in the system. 

To guarantee the compatibility with the constraints 
given by Eq. ([S]), we introduce two Lagrange multipliers 
\x a and n b into the mean-field Lagrange density Eq. |T]), 

K = L + ^ a |Va| 2 + M fc |^| 2 + (^+^)(l^| 2 + l^| 2 ). (7) 

Here the real parameters \i a and \i b can be identified 
as the chemical potentials of the corresponding compo- 
nents. Owning to Eq. ^ , we can conveniently introduce 
new notations: = ^i/^n a + n b , X lj = A^(n n +n b )/h, 
a = a'^TWjh, A, = A[(n a + n b ) 2 / 3 /h, Q = tl'/h, 
£ = £'/H, and A = A'/h. Reexpressing equation ([7]) in 
terms of <^>'s and substituting it into the Euler-Lagrange 
equation, we obtain a set of equations 

^ = E A -l^| 2 ^+^a|0 Q | 4/3 a +«^0 m -/i a Q , 

i ~dt = 5Z Abl l^| 2 ^ b + A b\M 4/3< l ) b + a(j>* a <j) m - fJ, b <j) b , 
= ^2\ m i\<j>i\ 2 <j) m + a<j) a (j) b -Q<i) g + £<j) m 
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= Y, X 3Mi\%-^m + ^g 



(8) 



where a phenomenological parameter 7 is introduced to 
characterize the decay of quasibound molecules. In terms 
of </>j's, the conservation relations (J6j) turn to be 



i<y 2 + i0j 2 + i</>/ = (i + £)A 



(9) 



(6) 



i&r + i0j^+i0„r = (i-*)/2, 

with S = (n a — n b )/{n a + n b ) characterizing the pop- 
ulation imbalance between fermionic atoms in different 
states. 



B. The definition of fidelity 

Now we are in the position to introduce a proper def- 
inition of fidelity for our system. As our system is a 
four-component system of Bose-Fermi mixture which is 
related to a graded unitary group SU(2 | 2), we need to 
define the fid elity carefully as it must obey several ba- 
sic properties [Ty] . Equation §§§ actually provides us the 
normalization condition 

|0 Q | 2 + |0 6 | 2 + 2|0 m | 2 + 2|0 9 | 2 = 1, 

which can be expressed as the following form 

{4>\F*{4>)F{cj>) 1 4>)=1, (10) 
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where (4> | denotes ((f)*, cf)* m , <jf ). One might 
think of a naive expression for the F'-matrix, F — 
diag(l, 1, \/2, V%)- However, because the relation of 
Eq. pop should be invariant under the transformation 
given in Eq. ([2]), just like that the conventional inner 
product in quantum mechanics is invariant under the 
U(l) transformation, the simplest correct expression of 
the F-matrix ought to be 



/ 


4>b 

\M 








\ 



















\<f>a\ 














V2 





V 











y/2 ) 



As a result, a nature definition of fidelity of a state la- 
belled by cf> with that labelled by <j)' is given by 

f(4>A')=\ (0|^(0W)I0'> I, (12) 

where the F-matrix was given in Eq. pip . Clearly, such a 
definition fulfills f(<P,U($ a ,$ b )q>') = /(f/^.^lW) = 
f((f>, 4>'), which means the phase transformation given in 
Eq. ([2|) does not vary the magnitude of fidelity; and the 
fidelity of a state with itself is always the unit /(</>, (j>) = 1 
which is just the normalization condition . 

III. COHERENT POPULATION TRAPPING 
STATES 

Now we consider the stationary states where we neglect 
the decay of quasibound molecules (i.e., 7 = 0). We 



know that the existence of stationary solutions of Eq. ( [8]) 
requires the system satisfies the adiabatic approximation 
condition. Once the adiabatic approximation is valid, 
i.e., difii/dt 0, the time-evolution equations ([5]) become 
algebraic ones for 0's, namely, 

Va4>a = A «l^| 2 0a + A a \<j) a \ i/3 <j> a + <*([)* (f> m , 

Vb<Pb = X] A 6il^| 2 06 + A b\<t>b\ A/3 <l > b + a##m, 
i^b 

(Va + Mfe)0m = X A "«l<^| 2 ^m + OL(j> a <t> b - +£<t>rn, 
i 

(/i a + /i b )0 fl = X A ^| 2< ^-^™ + A ^r ( 13 ) 

i 

Although it is difficult to find the exact solutions of the 
above equations, one can easily obtain a set of steady 
state solutions for Eq. P"3"|) by taking (j> m = 0. Such a 
state is called coherent population trapping (CPT) state 
which yields, 



_ (6-fr) + yJ(6-W)* + 2W(l + 8) 

l*2l a = m 2 -s, Kl 2 = ^-l0T, 

M a = Kb\<i>l\ 2 + K g \4>l\ 2 + A a \<&\^, 

H = Xabl^ + Xb^ + A^ 3 , (14) 

where = fi/a. The resonance condition corresponding 
to this solution is 



A = (A ab - \ ag Mf + (A afc - A bs )|^| 2 + (X ag + X bg \ gg M°\ 2 + A a \cjP a \^ + A|^| 4/3 . (15) 

I 



We consider a laser pulse with f2 — * 00 for t — ► and 
Cl — > for t — ► 00. For such a laser field, at the initial 
time (t = 0), we have 1 0° 1 2 = (1 + S)/2, |0°| 2 = 1-5/2 
and \<fP g 1 2 = 0, which implies there arc no molecules in 
the system at the initial time. At the final time (t — > 
00), |0°| 2 -> 5, |0 b | 2 -> and |0°| 2 - (1 - S)/2 for 
5 > 0; whereas |0°| 2 -► 0, \4> b \ 2 -> -8 and |0°| 2 -> 
(1 + (5)/2 for 5 < 0. With the help of initial values of 
\<f>® | 2 and their asymptotic values at final time, it is easy 
to find that those fcrmionic atoms, in the presence of 
their counterparts, can be converted into molecules if the 
CPT states can be followed adiabatically. The residual 
atoms can not be converted into molecules due to the lack 
of counterpart atoms. After the numerical calculation in 
next section, we will go back to study whether the CPT 
state can be followed adiabatically with the help of the 



useful concept, adiabatic fidelity. 



IV. ON CONVERSION EFFICIENCIES 

We know that there are two sorts of fermionic atoms 
6 Li and 40 K in the group of alkali-metal atoms. They 
have been converted into molecules in experiments suc- 
cessfully through the bare Feshbach resonance @, H, H, 
Hoj |. In those experiments, 60% to 80% of 40 K atoms and 
no more than 85% of 6 Li atoms can be converted into 
molecules. For 6 Li atoms, the atom-to-molecule conver- 
sion efficiency via bare Feshbach resonance is lower than 
that via STIRAP technique aided by Feshbach resonance, 
which is in contrast to the case for 40 K atoms. Due to 
the difference in atomic properties, the atom-to-molecule 
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conversion efficiency differs for different atoms even if the 
same technique is applied. 

Now we evaluate the atom-to-molecule conversion effi- 
ciency, respectively, for 6 Li and K atoms with concrete 
magnetic and laser fields. For these two sorts of atoms, 
we adopt the same time-dependent Rabi frequency, 



is very small. This is due to that \(j) m \ 2 is close to zero 
at any time in contrast to the case for 6 Li atoms. 



1 — tanh 



t-tr 



(16) 



where the parameters f2 max , an d t arc determined 
by the applied laser field that couples the two molec- 
ular states. In the numerical calculation, the detuning 
strength A is given by Eq. (fT5|) . We assume there are no 
molecules in the system at the initial time i.e., <fim.g = 
at t = 0. 




500 



FIG. 1: (color online) The time dependence of the population 
of particles for 40 K system for different 7. The parameters are 
5 = 0, X ab = 0.24, A a = A b = 0.16, £ = -4.4, fW = 200, 
to — 120, r = 40, and the other parameters are zero. Time 
is in unit of l/a and all other coefficients are in units of a, 
where a = 16.6 x 10" 29 J. 



For 40 K atoms, we know that the Feshbach resonance 
occurs at a magnetic field strength of 202. 1G and the 
resonance width is about 7.8G [l4j. Then we can get the 
atom-to-molecule coup ling strength a' = 16.6 x 10~ 39 J 
according to Ref. Here we choose that the mag- 

netic field is 201. 7G and particle density n a + rib is about 
10 20 m -3 . It is easy to obtain A a = 0.16a, A Q & = 0.24a 
and 8 — —4.4a. The time evolution of the correspond- 
ing population can be obtained by solving Eq. ((HI). We 
plot the numerical results in Fig. [TJ From this figure, we 



can see that the conversion efficiency 2|</> ff (t = 00) | 2 for 
40 K atoms is less than 60%, lower than that via the tech- 
nique of bare Feshbach resonance [l(| ■ The low conver- 
sion efficiency implies that the CPT state can not be fol- 
lowed adiabatically, which will be confirmed confidently 
by evaluating the adiabatic fidelity in next section. Com- 
paring the two panels in FigfTJ we can know that the in- 
fluence of 7 on the atom-to-molecule conversion efficiency 




FIG. 2: (color online) The time dependence of the population 
of particles for 6 Li system for different 7. The parameters 
are 5 = 0, \ ab = 0.0027, A a = A b = 0.055, £ = -0.1125, 
f^max = 200, to = 120, r = 40, and the other parameters are 
zero. Time is in unit of l/a and all other coefficients are in 
units of a, where a = 3.29 x 10~ 27 J. 



For 6 Li atoms, the Feshbach resonance occurs at two 
distinct strengths of magnetic field (the called narrow and 
broaden Feshbach resonance respectively). In our calcu- 
lation, we focus on the narrow Feshbach resonance for 
which the atom-to-molecule coupling strength a' is about 
3.29xlO _37 J. If the particle density n a + n b is about 
10 20 m~ 3 and the magnetic field is about 543. 6G, one can 
get A a = 0.055a, \' ab = 0.0027a and 8 = -0.1125a. 
Figure^] shows the time evolution of particle populations 
for the conversion of 6 Li atoms into molecules. From 
the top panel, we find that almost all of 6 Li atoms can 
be converted into molecules. For 6 Li atoms, the atom-to- 
molecule conversion efficiency via the STIRAP technique 
aided by Feshbach resonance is higher than that via bare 
Feshbach resonance @, Q , which is in contrast to the case 
for 40 K atoms. From the bottom panel, we can see that 
the atom-to-molecule conversion efficiency decreases dis- 
tinctly due to the existence of the quasibound molecular 
decay. 

The dependence of the conversion efficiency on the de- 
cay rate 7 of quasibound molecules for both 40 K and 6 Li 
is plotted in Fig. [31 Clearly, the influence of the decay on 
the conversion efficiency for 6 Li atoms is more distinct 
than that for K atoms, which confirms the previous 
interpretation. 



V. ADIABATIC FIDELITY FOR CPT STATES 

From the above numerical results, we can find that 
the atom-to-molecule conversion efficiencies of STIRAP 
technique aided by Feshbach resonance for 40 K and 6 Li 
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FIG. 3: The atom-to-molecule conversion efficiency versus the 
decay rate of the quansibound molecules for 40 K (left panel) 
and for 6 Li (right panel) system. The parameter choice for 
left (right) panel is the same as Fig. [1] (Fig. [2} . 



arc distinctly different although the CPT states are as- 
sumed to exist for the two systems. As we know that 
most of the atoms can be converted into molecules only 
when the CPT state is followed adiabatically, and the ex- 
istence of the CPT state does not guarantee the state can 
be followed adiabatically [HI]. The adiabatic properties 
for the atom-to-molecule conversion have been studied in 
Ref. Whereas, we can not apply this method to our 
system since all the possible solutions of the stationary 
equation |5]| must be known in the approach Q and it is 
difficult to obtain the other solutions beyond the CPT- 
state solution in our system. 
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FIG. 4: The time dependence of adiabatic fidelity for CPT 
state: (a) for 40 K, the same parameter choice as in Fig[T] (b) 
for 6 Li, the same parameter choice as in Fig[2] 

Recently, the fidelity was employed Q to character- 
ize the adiabatic condition for systems with atom-to- 
molecule conversion. The key point is that the value 
of the fidelity should close to unity if the system can 
adiabatically evolve in the CPT state. Now we evalu- 
ate the adiabatic fidelity f((j) (t), <j){t)) for our model, in 
which <p(t) is the wave function obtained by solving the 
dynamical Eq. ([5]) with 7 = 0, and (f>°(t) is the wave 
function corresponding to the CPT state which is given 



in Eq. (fT4f . The time evolution of the adiabatic fidelity 
for CPT state is plotted in Fig. Ufa) for 40 K system. One 
can see that the magnitude of fidelity is about 1 at the ini- 
tial time, which implies the system adiabatically evolves 
along CPT state. At the time 275/a, the fidelity dimin- 
ishes to the minimal value 0.35, which implies the system 
deviates away from the CPT state distinctly at that time. 
Although the fidelity begins to fluctuate later on, its fi- 
nal value is still no more than 0.75. The above analysis 
implies that the CPT state can not be followed adiabati- 
cally, so the atom-to-molecule conversion efficiency is not 
high for 40 K atoms. Fig.^b) shows the time evolution of 
the adiabatic fidelity for CPT state corresponding to the 
conversion of 6 Li atoms into molecules. Obviously, the 
fidelity for this system is very close to 1 at all the time, 
which implies that the system adiabatically evolves along 
the CPT state. This result is consistent with the fact 
that almost all of the 6 Li atoms can be converted into 
molecules. Since the fidelity for CPT state can well re- 
flect the atom-to-molecule conversion efficiency, one can 
improve the conversion efficiency by optimizing the pa- 
rameters of the system to achieve a higher adiabatic fi- 
delity. 



VI. SUMMARY 

With the help of mean-field Lagrangian density, we 
studied the conversion of two-species fermionic atoms 
into bosonic molecules via STIRAP technique aided by 
Fechbach resonance. We calculated conversion efficien- 
cies for 40 K and 6 Li systems respectively, and found 
that almost all of the 6 Li atoms can be converted into 
molecules, which implies that the STIRAP technique 
aided by Feshbach resonance is more effective than the 
bare Fechbach resonance for 6 Li atoms rather than 40 K 
atoms. We also compared the influence of the decay rate 
of quansibound molecules on the conversion efficiency for 
40 K and 6 Li systems respectively, and found that there 
is a big difference between them. The success of STF 
RAP technique not only requires the existence of the 
CPT state but also requires that the system can adia- 
batically evolve within such a state. 

The adiabatic fidelity was recently recognized to be 
a useful measurement for characterizing the adiabatic 
properties. Our analysis of the symmetry and the corre- 
sponding conservation law for the systems under consid- 
eration helped us to introduce an appropriate definition 
of adiabatic fidelity for CPT state. For 40 K system, the 
CPT state can not be adiabatically followed since the fi- 
delity was found to be less than 0.75 at final time, which 
is well consistent with the corresponding low conversion 
efficiency. In order to improve the conversion efficiency in 
40 K system, one should achieve a higher fidelity through 
optimizing the parameters of the system. Whereas, for 
6 Li system, the fidelity is very close to 1, hence the CPT 
state can be followed adiabatically, which is the reason for 
a high conversion efficiency. Our evaluation of the adia- 
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batic fidelity enable us to understand why the conversion Scholars and Innovative Research Team in University, 
efficiencies for 40 K and 6 Li are distinctly different. and NSFC Grant No. 10674117. 
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